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Abstract 


We analyse in detail the supersymmetric contributions to the decay of an 
extra Z boson in effective rank 5 models, including the important effect of 
D-terms on sfermion masses. The inclusion of supersymmetric decay chan¬ 
nels will reduce the Z' branching ratio to standard model particles resulting in 
lower Z' mass limits than those often quoted. In particular, the supersymmet¬ 
ric parameter space motivated by the recent Fermilab eeyy event and other 
suggestive evidence results in a branching fraction B(Z' —>■ e~^e~) ~ 2 — 4%. 
The expected cross sections and branching ratios could give a few events in 
the present data and we speculate on the connection to the three e'’'e“ events 
observed at Fermilab with large dielectron invariant mass. 
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I. INTRODUCTION 


Low energy supersymmetry provides the most successful solution to the naturalness 
problems which plague the standard model (SM). This has led to an extensive study of 
the experimental consequences of low energy supersymmetry in the literature and the next 
generation of colliders will explore a substantial fraction of the supersymmetric parameter 
space. If superstring theory is assumed to be the underlying fundamental theory at the 
Planck scale, then not only is it responsible for low energy supersymmetry, but the possibility 
exists that the SM gauge group may be larger at the TeV scale. This is because a consistent 
string theory requires large gauge groups (EgxEg or SO(32)). Depending on the particular 
compactihcation scenario, the resulting low energy N=1 supersymmetric theory often has a 
gauge group with rank larger than four. 

The simplest possible extension of the SM gauge group suggested by a gauge group of 
larger rank involves the introduction of one extra U(l) factor. This produces an extra neu¬ 
tral gauge boson, Z', in the particle spectrum. The low energy phenomenology of Z' bosons 
has been extensively discussed in the literature (see |]I| and references therein). Recently, 
particularly strong motivation for having the mass of the Z' below a TeV has been empha¬ 
sized by Cvetic and Langacker p[. However, most analyses of Z' physics do not discuss 
supersymmetric contributions in detail. In this work we will examine the supersymmetric 
decay channels of the Z' boson, including decays to neutralinos, charginos and supersym¬ 
metric Higgs bosons which are normally neglected. Our analysis will also include D-term 
corrections to the scalar masses which can have appreciable effects and are also not normally 
included Q 

The supersymmetric particle spectrum will be calculated by choosing typical values for 
soft masses and parameters which are consistent with recent suggestive evidence for 
supersymmetry, such as the eeyy event at the Fermilab Tevatron [^]. We prefer to focus 
on this region of parameters rather than survey the entire parameter space. This will 
lead to representative values of the Z' branching fractions in the presence of low energy 
supersymmetry. If a Z' exists then our analysis will provide a more realistic context in 
which to experimentally search for these particles. 

The outline of this paper is as follows: We begin in Sec. H by systematically writing 
down all the supersymmetric Lagrangian terms relevant for Z' boson decay in effective rank 
5 models. Using these terms we will derive the decay widths including the phase space 
factors. A parametrisation of D-term contributions to the supersymmetric scalar masses in 
effective rank 5 models is also briefly discussed and we comment on possible radiative U(l)' 
symmetry breaking scenarios. In Sec. HI we present a generic set of branching fractions for 
various Z' models. Using these values we show that the Z' mass bounds from direct limits 
at Fermilab are lower than those normally quoted. Final comments and the conclusion will 
be presented in Sec. IV. An Appendix also follows which summarises the mass contributions 


^While we were writing this paper some analysis of D-term contributions to scalar masses in Eq 
models was reported by E. Ma |^. The possibility of large mass shifts from the effect of Z' D-terms 
on sfermion masses has also been emphasized by J.Lykken [^. 
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for sfermions and Higgs bosons. 


II. Z' DECAY MODES WITH LOW ENERGY SUPERSYMMETRY 


The low energy theory resulting from the breakdown of extended gauge groups can 
include additional gauge bosons in the spectrum. Since the grand unihed gauge group Eg 
has rank 6 and the SM gauge group has rank 4 there can be at most 2 additional gauge 
bosons. For simplicity we will consider an effective rank 5 low energy theory with only one 
additional gauge boson associated with an extra U(l) and parametrised by 

Z'(6*) = Z.0 cos^ — Z^^sin^, (1) 


where 0 is a mixing angle and 

(i) occurs when Eg —>-80(10)xU(l)p, 

(ii) Z^ occurs when SO(IO)—>-SU(5)xU(l);;^. 

The orthogonal combination to (|^) is assumed to occur at the intermediate or Planck scale. 
When Eg breaks directly to a rank 5 group (SM x U(l)^) as in superstring inspired models 
via Wilson line breaking 0, the extra Z boson is denoted = y^S/S — y^S/S Z^. We 
will also consider the model with Q = rj — ti/2 which will be referred to as Z/. This model 
occurs when Eg breaks to SU(6)x 811(2)/ and is orthogonal to the p model. 

There are two types of mixing that can occur between the SM Z boson and an extra Z' 
boson. The hrst is associated with the fact that in general Z and Z' are not mass eigenstates. 
The mass squared matrix in the interaction basis has the form 


M| AM2 
AM2 M|, 


2gw g' E. T^iQm" 2 g'^ 


( 2 ) 


where gi,g 2 ,g' are the U(1)y, SU(2), U(1)' gauge couplings, gw = \jgl + gl and (0*) = 
Vi/\/2 is the vacuum expectation value (vev) of a Higgs held (pi with weak isospin Ta* 
and U(l)' charge Q'. This matrix can be diagonalised by an orthogonal matrix which is 
parametrised by a mixing angle p satisfying 


tan 20 


2AM2 
M| - M|, 


( 3 ) 


where Mz^ and are the mass eigenstates and 


Zi = cos p Z + sin p Z' (4) 

Z 2 = — sin 0 Z + cos 0 Z'. (5) 


Recent analysis of the LEP results indicate that this mixing must be small, p E 0.01, 
because the measured Z mass is quite close to the value predicted in the SM. Since the 
mixing angle in the effective rank 5 model is 

1^2 

0~2sin0iY(Q;cos2/3-Q'2sin20)^ ~ (6) 

Z2 Z2 
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TABLE I. 

U(l) charges of the fields 

in the 27 of Eq. 



Vs 

2^/TOQx 

2^/6 

2VT5Qn 

Q 

1/6 

-1 

1 

2 


-2/3 

-1 

1 

2 

d" 

1/3 

3 

1 

-1 

L 

-1/2 

3 

1 

-1 


1 

-1 

1 

2 


0 

-5 

1 

5 

H 

-1/2 

-2 

-2 

-1 


1/2 

2 

-2 

-4 

D 

-1/3 

2 

-2 

-4 


1/3 

-2 

-2 

-1 

S'^ 

0 

0 

4 
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where tan/? = v^jvx, the Z' boson needs to be well above the electroweak scale. Except for 
the decay mode Z' —> W^W~ this mixing will be nnimportant and for the most part we 
will ignore mixing effects and assume that Zi ~ Z and Z 2 — Z'. 

The other type of mixing is associated with the kinetic energy terms when there are two 
U(l) factors in the Lagrangian P-[Tl[|. This kinetic mixing will result in a shift of the U(l)' 
charges and depends on the matter content of the theory. It will also not be important for 
the results that we will obtain. 

We will assume that the matter superhelds reside in the fundamental 27 of Eg, which 
consists of the left handed helds 


27 = (Q, L, d\ v\ H, D\ H\ D, 


(7) 


where Q includes the left-handed quarks u and d and L includes the left-handed leptons 
u and e. The exotic matter superheld H is an electroweak doublet containing the exotic 
leptons N and E, while H'^ contains and W. The exotic superhelds D and are a pair 
of vector-like colour triplets, while is a standard model singlet. The SU(5) representations 
are 


10(g,n^e'=), 

5(L,dg, 

1 ( 3 ^). 

The hrst three representations of these make up the 16 of SO(IO), while the remainder form 
the 10 and 1. All the U(l) charge assignments of the helds contained in (J^) are given in Table 
|. In general the U(l)' charge of a held d* will be denoted g'($) = g^(d’) cos6* —g^($) sin 6*. 
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A. Fermion/sfermion sector 


The Lagrangian for Z' coupling to the fermions of the 27 is given by 

^ = 9h^ivf-af^5)fZ'^ (9) 

where / = ^ ^ ^ Dirac fermion and the vector, axial-vector couplings are 

«/ = i («'(/i) + Q'Ur)) = i l(0*(/i) + QMr)) - (Qx(Jl) + QAM) sin9] (W) 
“/ = ^ W'(/i) - Q'Ur)) = ^ [(Q+l/i) - 0 *(/r))cos9 - (QAIl) - 0x(/R))si“91 ■ (H) 

Note that in the above equations Q^fn) = —Qifl) we will assume that the U(l)' gauge 
coupling constant g'^ = {5/3)gf. The Z' decay rate into fermions is calculated to be 


./2 


r(z' - //) = 


■iP4 


-•>-*51 


1-4 


M|, 


1/2 


( 12 ) 


where Cf = 1(3) for leptons (quarks) |^. 

The Lagrangian describing the interaction of the Z' with sfermions is given by 

£ = ig'(vf ± af)fl j^jL,RZ'^ 


(13) 


where the +(—) is for left (right) handed sfermions respectively. The decay rate for Z' decay 
to sfermions is calculated to be 


./2 


9 

ml 

JL,R 


3/2 


r(Z' ^ flJL,R) = Cf^Mz> {vf ± I 1 - 4^ 


Ml 


(14) 


where the colour factor Cj is defined as in the fermion case and ^ is the Jl^r sfermion 
mass. In the case of non-negligible sfermion mixing (such as the top squark) we must work 
with the Lagrangian in the mass eigenstate basis 


C = i(J(v, ± afC0s2e/)fl^dJj^2Z'i, - ig'afSm2gf(f;dj2 + RdjAZ',, 


(16) 


where /i ,2 denotes the mass eigenstates (see Appendix). The decay width is similiar to (|^ 


with the appropriate couplings from (pAf ), except that the phase space factor in the case of 
different mass decay products is [1 — 2(m^ -|- m|)/M|, -|- (ml — m|)^/M|,]^A_ 


1. D-term co'ntributio'ns to scalar masses 

In general, breaking the U(l)' gauge symmetry with Higgs field vevs contributes to 
sfermion masses via U(l)' D-terms in the scalar potential |T2|-|^Q. Depending on the 
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number of Higgs fields which are used to spontaneously break the gauge symmetry the 
contribution to the scalar mass squared term has the form 

(16) 

i 

where the sum is over all Higgs helds which obtain vevs, (recall that {(jyi) = Vi/\/2). If in 
the pure rank 5 or p model the U(l)' symmetry is broken using only one SM U(l)' singlet, 
then the D-term contribution to the scalar mass squared may be written 

/2 

“ 5k|)(2VI5Q,(o)). (17) 

However in the effective rank 5 models which are parametrised by the angle 6, there are 
also orthogonal D-terms which come from breaking the extra U(l)" at an intermediate or 
Planck scale. Since these orthogonal D-terms terms can raise the sparticle mass spectrum to 
energy scales much greater than a TeV, we will assume that their contribution is negligible. 
This can be achieved by using a mirrorlike pair of U(l)'' charged SM singlets with charges of 
opposite sign to ensure D"-flatness. This amounts to assuming degeneracy of the mirrorlike 
Higgs soft masses at some high energy scale. Thus in the effective rank 5 models the D-term 
contribution to the scalar masses (assuming again that U(l)' is broken by one SM singlet) 
becomes 

+ Qs^’DQ'ia)- (18) 

In the special case of Eq breaking directly to a rank 5 gauge group we obtain the rj model 
and recover the result (|T7|) . The D-term contributions to the sfermion masses are listed in 
the Appendix. 

Note that the exotic superfields in the 27 will also contribute fermions and sfermions. 
However we will assume that these particles are heavier than the Z' mass and do not con¬ 
tribute to the Z' decay width. 


B. The Higgs sector 


The scalar Higgs doublets required for the spontaneous breakdown of the electroweak 
symmetry can be associated with the scalar doublet components of the superfields H and 
H^. In addition, the scalar component of the superheld is an electroweak singlet Higgs 
boson which is used to break the U(l)' gauge symmetry and give mass to the Z'. We will 
denote these Higgs bosons by 


$1 = 





(19) 


where ‘hi, 2,3 are the scalar components of H, and respectively. In the ground state 
the vacuum expectation values (vev) of the Higgs helds will be denoted by (0°) = Vi/\/2 
where i = 1, 2, 3. Note that only one generation of scalar doublets receive vevs. There are 
also scalar Higgs helds associated with the other two generations. These helds are referred 
to as “unHiggs” because discrete symmetries can be used to make sure that these other two 
generations of scalars do not acquire vevs i0- 
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1. Radiative U(l)' symmetry breaking 


While we can treat the vevs Vi as phenomenological parameters it is also interesting to ex¬ 
amine the consequences for effective rank 5 models if we assume that a radiative mechanism 
is responsible for breaking the U(l)' gauge symmetry. This is similiar to the analysis done 
by Cvetic and Langacker where it was shown that it is possible to achieve an Mz'/Mz 
hierarchy for general string models, without excessive fine tuning provided Mz' ^ ITeV 
(see also 0 ). We will consider the case of grand unification with effective rank 5 models 
parametrised by (|I[). If the radiative breaking of U(l)' is due to one SM singlet and the 
superpotential contains a term W = A<hi<I> 2 *h 3 , then at low energies we will assume that the 
scalar potential for the neutral Higgs bosons has the form 


V = + h.c) 

+ + Him? + 

+ + Q'MiP + W'lPP 


\4>' 


0 | 2\2 


( 20 ) 


where mi,m 2 ,m 3 are the soft supersymmetric Higgs masses and Ax is a soft trilinear pa¬ 
rameter. Given that V 2 = Vi tan /?, the minimum of the potential (pOD occurs at 


vi ~ 


vl ~ 


gl^il - taxA PY 

2ml 


-ml + ml taiY (3 


{Q'l + Q'2 tan^ (3) 


-mq 


g'V, 


{Ql + Q '2 iav? (3) 

« ’ 




( 21 ) 

( 22 ) 


where we have neglected the corrections from the potential terms involving A. The matrix 
elements of the Z — Z' mass mixing matrix become 


M|^ 

M|,^ 


(1 -|- tan^ /3) 


-4 


(1 — tan^ PY 

9 ' {Q'l - Q '2 Y) 

gw {l-iav?YY 


—ml -|- ml iaiY Y — 


(Q'l + Q'2 tan^ Y) 


Q'3 


m. 


-ml -|- ml tan^ Y 


(Q'l + Q'2tan^/9) 


mq 


2 m 3 -|- 8 ^ 


/2 


gl^ (1 - tan^ YY 


Q't(l-^) + QAl-^)tati"/3 

Vl V2 


X 


—mi -I- m 2 tan^ Y ~ 


(Q'l + Q '2 tan^ Y) 


mq 


(23) 

(24) 


(25) 


If tau/d = 0 and mf < 0 then the above expressions agree with those in [Q. In order to 
achieve a reasonable hierarchy between Mz and Mz' for negligible Z — Z' mixing </> we must 
have vj -C m^. If minima exist for ml < 0 in (^), then this can only occur when the charge 
combination Q'l -|- Q 2 tan ^/3 has the same relative sign as Q' 3 . If tan/3 ^ 1 this will happen 
when — 7 r /2 < 6 Y, —tt/S. For example, an exact numerical determination of the minimum at 
9 = —1.161 for the parameters tan/3 = 1.5, mi = m 2 = 100 GeV, m 3 = 500 GeV, Ax = I TeV 
and A = 0.02 yields an acceptable Z — Z' hierarchy with Mz' — 700 GeV and Y = —0.0057. 
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At values of 0 ^ —vr/S an extreme fine tuning of the soft breaking parameters is needed for 
a radiative mechanism to work. However if m\ > 0 in the potential (^) then we do not 
necessarily need the charge combination Qi + Q '2 f3 to have the same relative sign as 
Q 3 . In this case no hne tuning of the soft parameters is needed as v'l can be made small by 
cancellation of the positive terms against —m\ for all values of 9. Particular scenarios for 
achieving these various symmetry breaking potentials at low energies requires a complete 
renormalisation group analysis. 

Note also that V 3 ^ 00 when the charge Qg —>• 0, while Mz' — g'Q'sVs — The 

condition Qg = 0 actually occurs when 6 = 71 / 2 and signihes that there is no minimum to 
the assumed form of the potential (pOD, since the stabilising quartic term is proportional to 
Qg^ (of course, at the value 6 = 7r/2, ^g can no longer break the U(l)' gauge symmetry). 
Clearly there is a limit to how big one can tolerate V 3 in the above scenario and this will 
have consequences for the D-term contributions to the squark and slepton masses (see next 
section). 


2. Higgs decay modes 


After symmetry breaking the physical and unphysical charged Higgs bosons are given by 

= sin (3(l)f + cos (26) 

= cos — sin /30f (27) 


where tan/? = V 2 /V 1 and is the charged Nambu-Goldstone boson which is absorbed by 
the gauge boson. In the unitary gauge (G^ = 0) we will then obtain (pf = sin and 
02 = cos fdH^. The physical charged Higgs boson masses are obtained by diagonalising the 
mass mixing matrix given in the Appendix [|^, with the result 


m|± 


2XAv3 

-^ + 

sin 2(3 





(28) 


where A is the soft supersymmetry breaking parameter in the Higgs potential. The La- 
grangian for the Z' coupling to the charged Higgs states is 


c = ig\Q[ sin^ (3 - Q '2 cos^ l3)H+d>^H- Z'^ 

+ 9\Q'i + Q' 2 ) Axip cos(3MwiW^^^H- + W^^-H+)Z'^ (29) 


where Q\ = Q'{H) and = Q'{HQ. This gives rise to the decay widths 

>2 / \ 3/2 


( 2 \ ‘ 
l-4^j 


487r' 

9 '^ 


(30) 


r(Z' ^ W^H^) = ^^(Q'l + Q'^Y siiA (3 cos^ (3Mz. 


1927r 




Ml 


Ml 


X 


\ 


^ _ ^ (M^ + m^±) ^ {Ml, - ml±Y 


Ml 


Ml 


(31) 















where is the mass of the Higgs boson. This agrees with the result in the limit 
mH± -C Mz'. 

While the Z' boson has no direct couplings to the W-boson, the small mixing angle 
discussed earlier in Eq.(^ induces a Z'W~^W~ coupling via the small Z admixture. The 
decay amplitude into W~^W~ pairs is calculated to be 


r(Z2 ^ W+W-) = 


1927r 


cos^ 9w sin^ (j)Mz. 




1-4^ 

. Ml) 




Ml 


Ml. 


(32) 


Recall that in the limit ^3 3> uig, the mixing angle (H) is 0 oc and so the width 

(^) in the equivalence theorem limit ^/s 3> Mw becomes 


/2 


T{Z' W+W) ~ T{Z' G+G-) = ^(Q'l cos^ p - Q'^ sin^ pfMz'. 

4&7T 


(33) 


Note that in ( P^) the Z — Z' mixing angle 0 cancels the usual enhancement factor, 

that one normally expects in the equivalence theorem limit. Similarly, in the limit 3> ^ 1,2 

the equivalence theorem dictates that M] 


/2 


T{Z' W^H^) ~ T{Z' G^H^) = Pcos^ pMz'. 

487r 

In the neutral Higgs boson sector, the Higgs fields are written as 


(34) 




(35) 


where pm, pu are the real and imaginary parts of 0°. Of the six neutral Higgs degrees of 
freedom, two (G°, G'°) will act as the Nambu-Goldstone bosons and will be absorbed by the 
Z and Z'. The remaining four degrees of freedom will become the physical neutral Higgs 
bosons which consists of three scalars, 2 3 and one pseudoscalar P°. In the unitary 
gauge (G, G' = 0) we will obtain PD|J2T| ,P| 


(.0 = !^po 
AT 

0^2 0 
N 

= E 

i=i 




(36) 

(37) 

(38) 

(39) 


where N = \Jv 1 V 2 + nfni + ^ 2^3 and U is the inverse of the matrix that diagonalises the 
scalar Higgs {H^) mass term. The Lagrangian relevant for Z' boson decay to neutral Higgs 
bosons is 
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( 40 ) 


£ = 2g'Mz E(Q'i cos^Uu - Qfs sin/3 U 2 ,)Z'^Z‘‘H^ 

i=l 

3 

“ 9 ^^{vsQ'isin (3 Uu + V 3 Q 2 COS U 2 i + vQ'^sin f3 cos f3U3i)Z'^H°d^Po, 


(41) 


2 = 1 


where v = \fvi + V 2 - The Z' —> ZHf decay width may be obtained in a similar manner to 
(^) and is given by 


/2 


r(Z' ^ ZH^) = ^{Q[ cospUu - Q'2smf3U2ifMz' 

dSvr 


(5M| - mjro) (M| - m^o)^ 
1 + 2-—^^ + - 


Ml 


Ml 


X 


\ 


(M| + m^o) (M| - m^o)2 
1 - 2 -^ + - 


Ml 


Ml 


(42) 


where ttiho is the mass of the scalar Higgs boson In general the 3x3 scalar Higgs mass 
mixing matrix is complicated (see Appendix) and mnst be nnmerically evalnated. However, 
in the limit that ^3 ^ ni ^2 the state 0 ° deconples and the only mixing occnrs between 0 ^ 
and 02 - If this mixing is parametrised as 




cos a — sm CK 
sin a cos a 




h ^2 r 


(43) 


then onr resnlts agree with those in the literatnre |r8|,p]|. 

Similarly, the Z' decay into scalar-psendoscalar Higgs bosons is 


r(z' ^ hf°Po) = 


g'^ 

4^i^ 

X 1-2 


{v^Q'i sin (5Uii + V 3 Q 2 cos (5 U 2 i + nQs sin 0 cos 0 U^ifMz 
{m]jo+mil) 


Ml 


+ 




(44) 


where mp^ is the mass of the psendoscalar Higgs boson (Po)- The mass of Pq is mnch easier 
to obtain in the exact limit. Using the mass matrix in the Appendix one hnds that 



2XAv3 
sin 20 


1 + 


4n| 


sin^ 20 


(45) 


If one again considers the limit ns S> ni ^2 and the mixing (^31), then the decay widths for 
Z' P°P °2 agree with the expressions obtained in 


C. The neutralino sector 

In the Eq models that we are considering there will be two extra nentralinos in addition to 
the fonr fonnd in the minimal snpersymmetric standard model (MSSM). An extra nentral 
gangino comes from the extra Z' boson and an extra nentral Higgsino comes from the 
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Higgs boson which is used to break the U(l)' gauge symmetry. In the interaction basis 
{-iB, -iHds, -iB\ $ 2 , $ 3 ), the 6 X 6 neutralino mass matrix is given by 


Ml 

0 

0 

-^givi 

\glv 2 

0 


0 

M 2 

0 

lg2Vi 

-\g 2 v 2 

0 


0 

0 

M' 

Q'lg'vi 

QWv2 

QW'^s 

(46) 

-^ivi 

\ 92 VI 

Q'lg'vi 

0 



\g 1 V 2 

-\g 2 v 2 

Q2g'v2 


0 



0 

0 




0 



where Mi, M 2 , and M' are the soft gaugino masses for B,Wz and B', respectively. The 
supersymmetric Higgsino mass terms come from the trilinear superpotential interaction 
y<|)i$ 2<|)3 which gives rise to a /i term with /r = Xvs/\/2. Q The neutralino mass eigenstates 
(A^i, N 2 , ■ ■ ■, Nq) are obtained by diagonalising the mass matrix (^ 6 [) . A general analytic for¬ 
mula for the neutralino masses does not exist. However, in certain limits such as ^3 3> Ui, V 2 
approximate expressions can be obtained p3. Of course, one can simply diagonalise the 


mass matrix (|46|) numerically, and this will be done for the results in the next section. 

In terms of the mass eigenstates {Ni, N 2 , - ■ ■, Nq) we may parametrise the coupling be¬ 
tween the neutralinos and the Z' boson as 


£ = JlaNirtAK 


(47) 


where the coupling constants gij are obtained from the diagonalisation of the neutralino 
mass matrix (^). 

The expression for the Z' decay width into neutralinos can be obtained using the La- 
grangian ( 0 )- The decay widths are 


V{Z' ^ NiN,) = ^Mz' 




127r 


1 - 


{mf + m3) {mf - 


2M|, 


2M|, 


^mimj 

Ml 


X 


\ 


(jUi + mj)‘^' 

Ml 


[mi — mj)‘^ 


Ml 


(48) 


where the m* refer to the neutralino masses. This agrees with the expression in |^. The 


coefficients gij will be determined numerically when we calculate the neutralino branching 
fractions in the next section. 


D. Chargino sector 

Since the Z' boson and the singlet Higgs supermultiplets are electromagnetically 
neutral they do not contribute any extra particles to the chargino spectrum. Consequently 
the chargino mass matrix remains the same as in the MSSM, namely 


^One can assume that the broken U(l)^ gauge symmetry generates an effective g term as advocated 
by Cvetic and Langacker Q. 
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( 49 ) 


M 2 \/2My[/ sin P 
\/2M]y cos /3 —/i 


except that /i = \v^/\/2 as defined above. The Z' can of course couple to the charged 
Higgsinos which leads to the chargino Lagrangian 


^ = \9' Y. + aij'j5)CjZ'^ ( 50 ) 

^ i,j=l 

where Ci are the chargino mass eigenstates and 

= Q[ sin^ (j)_ — Q'2 sin^ ( 51 ) 

flu = Q'l sin^ 0 _ + Q'2 sin^ ( 52 ) 

1^12 = 1^21 = Qi sin 0 _ cos 0 _ — 6Q'2 sin 0 + cos 0 +, ( 53 ) 

ai 2 = 021 = Q[ sin 0 _ cos 0 _ + 5(^2 sin 0 + cos 0 +, ( 54 ) 

V22 = Q'l cos^ 0 _ - Q'2 cos^ 0 +, ( 55 ) 

022 = Q'l COS^ 0 _ + Q'2 COS^ 0 +. ( 56 ) 


In (il-dl), (j)± are the angles of the unitary transformation matrices used to diagonalise 
the chargino mass matrix m and 6 = nin 2 , where ni = sgn(m(^±) and n 2 = sgn{mQ±). 
The Z' decay rate into chargino pairs is calculated from the Lagrangian (|5I1|) to be 


r(z' ^ c,Cj) = 


-J2 


dSvr 


-M71 


12 , 2 i(, (m^ + mj) 


2Ml, 


2M% 


I Q/ 2 2 1 


X 


\ 


1 - 


{rrii + nijY‘‘ 

Ml 


1 - 


(m* — rrijY 

Ml 


( 57 ) 


In the limit that Mz' 3> m^. the decay width (^) reduces to a form similar to the result 
(O) that was obtained earlier in the fermion sector. 


III. Z' PRODUCTION CROSS SECTION AND BRANCHING FRACTIONS 


At the pp Tevatron collider the Z' production mechanism is due to the Drell-Yan process. 
If the u and d quark couplings to the Z' are obtained from (^, then the expression for the 


Drell-Yan production differential cross section is p4| 


da , 
dy 


{AB Z'X) = Y i^q + «?) q{xa) q{xb) 


/2 


( 58 ) 


where y is the rapidity of the Z' in the AB c.m. frame, k ~ 1 -|- ^as(M|,) is a QCD 
correction factor, Xa{xb) is the momentum fraction of q{q) in A{B) and q{x) is the quark 
distribution function. Similarly, the Z' production cross section at the planned 14 TeV pp 
LHC collider is also obtained from (jSSf) . In Fig. |I|the cross section for Z' production in the 
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M^, (GeV) 


FIG. 1. (a) The Z' production cross section for various effective rank 5 models at the Fermilab 
Tevatron with ^/s = 1.8 TeV. (b) Similarly for the LHC with ^/s = 14 TeV. 




TABLE II. Branching fractions of all possible Z' decay channels in a supersymmetric framework 


for various Z' models. 


Z' decay channel 

Zi 

z^ 

Z-q 

e'’“e“ 

0.0391 

0.0280 

0.0171 

/i+/X“ 

0.0391 

0.0280 

0.0171 

+ 

1 

0.0391 

0.0280 

0.0171 

VeVe. 

0.0782 

0.0280 

0.0890 


0.0782 

0.0280 

0.0890 


0.0782 

0.0280 

0.0890 

uu 

0.0000 

0.0839 

0.0820 

cc 

0.0000 

0.0839 

0.0820 

it 

0.0000 

0.0553 

0.0540 

dd 

0.1174 

0.0839 

0.0513 

ss 

0.1174 

0.0839 

0.0513 

bb 

0.1174 

0.0839 

0.0513 

J2UiUj 

0.0000 

0.0000 

0.0000 

Ed*dj 

0.1746 

0.0000 

0.0000 

Y.e*ej 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.1280 

H+H- 

0.0048 

0.0021 

0.0003 

EPoH! 

0.0110 

0.0038 

0.0010 


0.0018 

0.0102 

0.0039 

EZHf 

0.0075 

0.0500 

0.0206 

W+W- 

0.0028 

0.0062 

0.0155 

ECiCj 

0.0182 

0.0759 

0.0329 

ENiNj 

0.0753 

0.2090 

0.1077 


'll), X and r] models at both the Tevatron and the LHC are depicted. The slight variation in 
the cross section for each Z' model in Fig. |I] is dne to the 9 dependence of the qnark-gange 
boson conplings ( ([lOD and (|Tl|)). 

Direct limits on pp ^ Z' ^ e~^e~ at the Tevatron can place a lower bonnd on the Z' 
mass depending on the valne of the branching ratio B{Z' —>■ |^. In order to estimate 

the e^e“ branching ratio in a snpersymmetric framework, all possible decay channels (snch 
as decays to sparticles) need to be inclnded. While there are many nnknown parameters 
in the snpersymmetric standard model it is possible to constrain parameters by reqniring 
consistency with other indirect stndies of snpersymmetric parameter space. For example, in 
the recent snpersymmetric interpretation of the eeyy event, the snpersymmetric parameters 
/i, tan/?. Ml and M 2 are constrained 0- The parameter ranges of these variables will in 
tnrn constrain the nentralino and chargino mass spectrnm. 

Similarly, by assnming typical ranges of snpersymmetric parameters which are reqnired, 
for example, in the constrained minimal snpersymmetric standard model the branching 
fractions of all possible decay modes can be estimated. In Table we list the branching 
fraction of all possible Z' decay channels in various Z{6) models for the following choice of 
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supersymmetric parameters: 


tan (3 = 1.5, /i = —50 GeV 
Ml = 80 GeV, Ma = 100 GeV, M' = 300 GeV 
m = 500 GeV, A = 500 GeV 

where we have assumed common scalar mass m and A terms. In addition we also £x the Z' 
mass to be Mz' — 700 GeV, which sets the scale for all the kinematically allowed decays. 
Thus, in Table when the branching fraction B(Z' —> f*fj) = 0 it is because the squark 
and slepton masses are too heavy to be kinematically allowed. Similarly only the sufficiently 
light Higgs bosons, neutralinos and charginos are included. All exotic particles such as the 
vectorlike quarks are also assumed to be very heavy. 

In Fig. ^ we show the squark and slepton masses including the D-term contributions. 
The general effect of the D-term contributions is to make the squarks and sleptons heavy 
since the scale of the D-terms is set by Mz' (or ^ 3 ). However, in Fig. ^ there are special 
values of the mixing parameter 6 where the D-term contributions cancel and the squark and 
slepton masses become light. These are the only points that could be consistent with squarks 
and sleptons that would have appeared already at Fermilab [^. Such sensitivities are very 
encouraging from the point of view of extracting information about new physics from limited 
data. Beyond these values the squared masses become negative and this signals the onset 
of charge and colour breaking minima. It may be possible that radiative corrections can 
stabilise the vacuum but the analysis of these corrections is beyond the scope of this paper. 
When the Z' mixing parameter 6 approaches ± 7 r /2 the squark and slepton masses become 
unacceptably large. This is because the D-term contribution to the scalar masses from ^3 
(Eq. (PIBD) is Am^ oc and when Q'^ —>• 0 we have Am^ —>• ±cx). Thus, in effective 

rank 5 models with a U(l)' symmetry breaking potential of the form (^Of) one can exclude 
y-like models (or regions near 6 ~ =t 7 r/ 2 ) because the squark and slepton masses become 
large and negative leading to charge and colour breaking minima (where we have neglected 
radiative corrections). 

As far as determining a lower Z' mass bound from collider experiments, the main conse¬ 
quence of accurately including all supersymmetric decay channels is that the branching frac¬ 
tion into leptons and quarks is reduced. This lowers the Z' mass bound from that normally 
quoted in the literature. The branching fraction from Table |T| is in the range 2 — 4%, 
where typically in non-supersymmetric models one obtains 3 — 6 % |^. If we use the direct 
limit from a recent combined analysis of GDF and DO data {a x B{Z' e~^e~) < 0.28pb) 
P5| , then for the branching fractions listed in Table we obtain Mz> ^ 360, 370, 360 GeV 
for Zj, Z^ and respectively. 

The importance of various decay channels can be determined from how the branching 
fractions in Table ^ vary as a function of the Z' mixing parameter 9. This dependence is 
shown in Fig. ^ for rh = 500 and 1000 GeV. We have excluded the regions corresponding to 
negative squark and slepton masses which are near 9 = ± 7 r/ 2 . As the soft mass parameter 
rh becomes larger the charge and colour breaking regions shrink. Near the excluded regions 
(which correspond to and Zj for m = 500 GeV) the Z' branching fractions to squark 
and sleptons are maximised because their masses become light enough to be kinematically 
allowed. As 0 —0 the D-term corrections make the squarks and sleptons kinematically 
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^ (TeV^) mg (TeV^) m| (TeV^) m~ (TeV^) 



e 

FIG. 2. The squark and slepton masses squared as a function of the Z' mixing parameter 6. 
Mixing is only important for the top squarks and the solid line in the up squark plot represents ti^ 2 - 
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Branching Fraction Branching Fraction 



I I W+W, W^H~, 
H^W~, 'LZHf 

IM H+H~, ZP°Hf 

□ v^, V, 

□ e, n, T 
I I u, c, t 
I I d, s, b 


FIG. 3. The Z' decay branching fractions as a function of the mixing parameter 6 for a repre¬ 
sentative set of supersymmetric parameters as defined in the text. The excluded regions correspond 
to values of the squark and slepton masses that lead to charge and colour breaking minima. The 
lower figure has m = 1 TeV 


17 








inaccessible. A significant neutralino and chargino branching fraction occurs in models 
which are ■^-like (or 0 ~ 0). Specifically, decays to the lightest neutralino can be as large as 
10% and this would contribute greatly to the Z' invisible width. As 6* —±7r/2 the branching 
fraction to neutralino and chargino pairs gets smaller while the neutrino branching fraction 
becomes larger. Decays to Higgs bosons become non-negligible for 0^0. One should 
also note that the branching fraction to leptons and quarks remains fairly constant. In 
particular the branching fraction to light quarks (u,c,d,s,b) is at most 0.45. Thus for light 
quarks a x B{Z' —> qq) ~ 0.25pb where we have assumed mz' — 700 GeV. This is too small 
to be observable in the inclusive jet cross section. Note that the branching fraction to up 
quarks vanishes for Z/-type models. This is a consequence of the fact that the up-quark Z' 
coupling becomes zero. 

It is also amusing to note that both CDF and DO report events with dielectron 
invariant mass ^ 500 GeV and CDF also has a second event with invariant mass ~ 350 
GeV. If one assumes that these high energy events are due to the decay of a Z' boson, then 
given our range of B(Z' —> e'''e“), this would correspond to a Z' mass ~ 600 — 700 
GeV. One would also expect these events to be backward, i.e. cos^* < 0 where d* is the 
angle of the outgoing e~ with respect to the quark in the qq centre of mass frame p^ , |3^ . In 
Fig. ^ we plot the forward-backward asymmetry, ApB as a function of the mixing angle 6. 
A very large asymmetry is expected for Zj type models. It is tantalising that the two events 
reported by CDF with invariant masses of 350 GeV and 504 GeV have cos 6* values of -0.14 


and -0.27 respectively 


If we require that the probability of observing two backward 


events be at least 1/e then one would need A 


FB 


< 


- 0 . 2 . 


IV. CONCLUSION 

We have examined in detail the supersymmetric contributions to the decay of an extra 
Z boson in effective rank 5 models parametrised by the angle 6, including the important 
effect of D-terms on sfermion masses. The supersymmetric particle spectrum was chosen so 
as to be consistent with other analyses such as the recent eeyy event and the constrained 
minimal supersymmetric standard model. The main effect of including these contributions 
is that it reduces the Z' branching fraction into lepton pairs. This, in turn reduces the lower 
bound on the Z' mass obtained from direct limits at the Tevatron. 

While the Z' decay to neutralino and chargino pairs is always non-negligible, this is not 
necessarily the case for squarks and sleptons. When the squarks and sleptons are dominated 
by positive D-term contributions their masses can become very large and consequently they 
are no longer a kinematically accessible Z' decay mode. However for special ranges in 
parameter space the D-term contributions become negative and give rise to light squark 
and slepton masses. At these values the branching fraction of squarks and sleptons is 
non-negligible. As 9 approaches ±7r/2 various squark and slepton masses squared become 
negative leading to charge and colour breaking minima. Thus depending on the size of soft 
squark and slepton mass parameters, effective rank 5 models with regions of 6 near ±7r/2 
produce unacceptably large D-term contributions. If one also assumes a U(l)' symmetry 
breaking potential (PO]) then it is possible to achieve a radiative breaking mechanism without 
any fine-tuning of the soft parameters. 
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IG. 4. The forward-backward asymmetry for e'^e pairs a 
function of the Z' mixing parameter 6 with Mz> = 700 Ge 




Finally, the forward-backward asymmetry of e’''e“ pairs becomes quite large at the Zj 
model and the two reported CDF events are consistent with models of this type. The 
measurement of Aps will provide an interesting test for Z' models at future colliders. 
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APPENDIX; 

This appendix summarises the D-term contributions to the sfermions and the neutral 
scalar Higgs boson masses. 


1. Sfermion masses 


The sfermion mass matrix is parametrised as 


Mj = 


Mh Mlp 

^LR ^RR 


(Al) 


Dehning the U(l)' D-term contribution to be m%, = + Q 2 V 2 + ^^e mass¬ 

mixing matrix elements for the up squarks ul,r are given by 


^LL — Af&z, + + (o 


-xw)Ml cos 2(3 + Q'u^rhy 


Mrr = Ml^ +ml+ ‘^xwMl cos 2(3 + 

m; 


u 2 
LR 


= mu{Au - ficotp), 


and for the down squarks dp^R 


Mfi = + (-o + oXw)Ml cos 2(3 + 


dL 


Mrr = + -xwMl cos 2(3 + 

T 2 

^LR = -/itan/3) 

where xw = sin^ 6w Similarly for the cr^r sleptons we obtain 


(A2) 

(A3) 

(A4) 

(AS) 

(A6) 

(A7) 
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^LL = +ml + + xw)Ml cos 2/3 + 

(A 8 ) 

+ml- xwMl cos 2/3 + 

(A 9 ) 

^LR = rne{Ae - /itan/ 3 ), 

(AlO) 

sleptons 


MIl = cos 2/3 + 

(All) 

M^rr =Ml + 

(A 12 ) 

Min = 0 . 

(A 13 ) 


The sfermion mass eigenstates are given by 


/i 

fi 


cos 9 j sin 6 »^- \ ( h 
-sinOj cos 9 j J [ 


(A14) 


This basis is really only important for the top squarks where the mixing term is non- 
negligible. 


2. Higgs masses 


V — ~ + h.c.) 

(AI6) 

+ A^(<hl<hi<h^<h2 + <hl<hi<h^<h3 + <h^<h2<h^<h3) + (y “ A^)|<hl<h2p 

(AI6) 

+ - 4>^4>2)2 + + Q'2H^2 + (A17) 

where we are assuming a superpotential term W = A<hi$ 2 ‘h 3 with defined as in ([I^) and 

we have written the D-terms for the more general effective rank 5 model. The neutral Higgs 
bosons iff masses directly receive D-term contributions from the spontaneous symmetry 
breakdown of the extra U(1)A The mass mixing matix is given by 

Mjjo — - 

Bivl + \Av2V^/vi B2V1V2 — XAv^ B3V1V3 — XAv2 

B2V1V2-XAV3 B4V2 + XAV1V3/V2 B5V2V3-XAV1 
B3V1V3-XAV2 B5V2V3-XAV1 Bev'^ + XAV1V2/V3 _ 

(AI8) 

where 

Bi = 2^^! A ^ 2 ) + "^Q'l 9'“^ 

(A 19 ) 


B2 = 2A^ — ^(S'l + 92) + ‘^Q'iQ'29'^ 

(A 20 ) 
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The Higgs boson masses are obtained from the Higgs potential Em 






B3 = 2X^ + 2 Q\Qy‘^ 

(A21) 

-®4 = g'i) + 2Q2 9 '^ 

(A22) 

H5 = 2A2 + 2Q'2Q39'^ 

(A 23 ) 

B, = 2 Q'^^g'‘^ 

(A 24 ) 


For completeness we also list the pseudoscalar and charged Higgs boson mass matrices |]l| 


M\,o 


\Av^ 

2 


Vi/V2 1 U 1 /U 3 

1 V2/V1 V2/V3 

V1/V3 V2/V3 V1V2/VI 


(A 25 ) 


{g 2/2 - X^)vl + XAV 1 V 3 /V 2 {gl /2 - X^)viV 2 + XAv^ 
{gl /2 - X^)viV 2 + XAv 3 {gl /2 - X^)vl + XAV 2 V 3 IV 1 

where we have used the definitions in Sec.II. 



(A 26 ) 
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